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ABSTRACT 



Aims. We study the environmental dependence of local luminous infrared galaxies (LIRGs) and ultraluminous infrared galaxies 
(ULIRGs) found in the Sloan Digital Sky Survey (SDSS) data. 

Methods. The LIRG and ULIRG samples are constructed by cross-correlating spectroscopic catalogs of galaxies of the SDSS Data 
Release 7 and the Infrared Astronomical Satellite Faint Source Catalog. We examine the effects of the large-scale background density 
(E5), galaxy clusters, and the nearest neighbor galaxy on the properties of infrared galaxies (IRGs). 

Results. We find that the fraction of LIRGs plus ULIRGs among IRGs (/(U)Lirgs) and the infrared luminosity (Lir) of IRGs strongly 
depend on the morphology of and the distance to the nearest neighbor galaxy: the probability for an IRG to be a (U)LIRG (/(U)lirgs) 
and its Lir both increase as it approaches a late-type galaxy, but decrease as it approaches an early-type galaxy (within half the virial 
radius of its neighbor). We find no dependence of /(U)lirgs on the background density (surface galaxy number density) at fixed stellar 
mass of galaxies. The dependence of /(U)lirgs on the distance to galaxy clusters is also found to be very weak, but in highest-density 
regions such as the center of galaxy clusters, few (U)LIRGs are found. 

Conclusions. These environmental dependence of LIRGs and ULIRGs and the evolution of star formation rate (SFR)-environment 
relation from high redshifts to low redshifts seem to support the idea that galaxy-galaxy interactions/merging play a critical role in 
triggering the star formation activity of LIRGs and ULIRGs. 

Key words, galaxies: active - galaxies:evolution - galaxies:formation - galaxies:interactions - galaxies :starburst - infrared: galaxies 



1. Introduction 

Since the launch of the Infrared Astronomical Satellite 
(iNeu^ebauer et all 11984 IRAS) in 1983, a great deal of ef- 
fort has been made to understand the nature of luminous in- 
frared galaxies (LIRGs; 10" < Lir < lO^^^o) and ul- 
traluminous infrared galaxies (ULIRGs; Lir > IO'^Lq); (1) 
what triggering mechanism is responsible for the enormous in- 
frared (IR) luminosity (internally or externally driven?); (2) 
. what determines the relative importance between star forma- 
' tion (SF) and active galactic nuc l ei (AGN) activity within 
■ (U)LIRGs (see Sanders & Mirabel' '1996; 'Genze l & Cesarskvl 
|2000; Lonsdale et al. 2006; Soifer et al. 2008 and references 
' therein). For example, most ULIRGs are found to be interacting 
, systems between two or mo re late-type galaxies (e.g. jKim et alJ 
I2OO2L IVeilleuxetal]l2002h . Many LIRGs are also found asso- 
ciated to interacting systems, but a significant fraction of them 
do not show defi nite featur es of major me r gers (J^id.a .2()04; 
Wang et alJ 120061 see also IShi et alJ l2006l l2009nElbaz etalJ 



and extensively studied (e.g ., Lewis et al. [2002 : Gomez et al 



'2003; Kauff^mann et al. 2004; Balogh et al. 2004; Tanaka et al 
2004; Owers et al , 2007; Park et a l. 2007; LietaL 200| 
iPark&Hwand 120 09; Hwang &^3|200: 



Par k & Choi - .. _ 

Blanton & Moustakasl2009l;lLee et an2010 ). However, there are 



few studies focusing on how (U)LIRGs, showing much larger 
SFA compared to normal galaxies, are aflrected by the environ- 
ment. 



200f iLotz et al. I I2OO8I: lldeue et alJl2009l) . Thus, minor mergers 



and galaxy-galaxy interactions may also play an important role 
in triggering the star formation activity (SFA) of LIRGs. 

As the Sloan Digital Sky Survey (SDSS; lYork et alJ 

I2OOOI) and Two Degree Fi eld Galaxy Redshift Survey (2dF- 
GRS; ICoUess et al.l l20()Th have produced unprecedentedly 
large photometric and spectroscopic data of nearby galax- 
ies, the important role of the environment in determining 
the SFA of galaxies in the local universe has been revealed 



iGotd (l2005h investigated the optical properties of IR galax- 
ies in the SDSS Data Release 3 (DR3) and found that more 
luminous IR gala xies tend t o be lo cated in lower density re- 
gions. Later, Zauderer et a D (I2007h studied the environment 
of local ULIRGs at z < 0.3 using the spatial correlation 
amplitude parameter Bgc. The environment of most ULIRGs 
was found to be similar to that of the field, but a small 
number of ULIRGs was found in high-density regions such 
as galaxy clusters. They also found no significant difference 
in the Bgc distribution of local ULIRGs with those of lo- 
cal Seyferts, local QSOs, and intermediate-z QSOs, Recently, 
Wang & Rowan-Robinson (2010) , usi n g the Imperial IRAS- 
FSC redshift catalog ( Wang & Rowan-Robinsonll2009l) . inves- 
tigated the spatial clustering of galaxies as a function of sev- 
eral galaxy properties. They found that cirrus-type galaxies (dust 
emission due to the interstellar radiation field) show a stronger 
clustering than M82-like starbursts (interaction-induced star- 
bursts), which indicates that the latter preferentially resides in 
low-density environments. They also found a correlation be- 
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tween the correlation strength and the IR luminosity for M82- 
like actively star-forming galaxies except in high-density regions 
such as the center of galaxy clusters. 

The environment of (U)LIRGs beyond the l o cal un iverse 
has recently started to be explored. iFarrah et"al] (l2006l) stud- 
ied a spatial clustering of ULIRGs at 1.5 < z < 3, and 
found that these ULIRGs clust er more stro ngly than other galax- 
ies at these epochs. Similarlv. iGilli et aTl (|2007) suggested that 
24/im-selected, star-forming galaxies at z ~ 1 in the Great 
Observatories Origins Deep Survey (GOODS) appear to be more 
clustered than optically selected galaxies at similar redshifts, and 
that more luminous IR galaxies have a larger correlation length. 
Recentlv. lMarcillac et al.l (l2008h investigated the environment of 
LIRGs and ULIRGs at 0.7 < z < 1 in the Extended Groth Strip 
(EGS). They found that the local environment of (U)LIRGs is in- 
termediate between those of red and blue galaxies. If they com- 
pare (U)LIRGs with non-IR galaxies whose stellar mass, color, 
or luminosity are similar to those of the (U)LIRGs, they see no 
significant difference in the environment between the two. 

Since LIRGs and U LIRGs domi nate the S FA at high red- 
shifts (e.g.. lLe FlocTet al. 2005: Magnelli et al. 2009) . the un- 
derstanding of the nature of (U)LIRGs is closely related to 
the study of the evolution of the star fo rmation rate (SFR)- 
environment relation. For example. lElbaz e t al. (2007), using the 
GOODS data, discovered that the SFR-density relation observed 
locally is reversed at z ~ 1 , in the sense that the "average" SFR of 
galaxies over a given volume increases with the local galaxy den- 
sity (see also 'Cooper et al."2008- 'Ideue e t al.l [2009HTran et all 
[2010). On the conti-ary, Patel et al. (2009) observed the galax- 
ies in a galaxy cluster (RX J0152.7-1357) at z = 0.834, and 
found that the SFR-density re lati on at that redsh ift is not similar 
to those in lElbaz et all ( l2007h and lCooper et 311 (2008), but simi- 
lar to that in the local universe (see also Feruglio et al. 2009). It is 
noted that the SFR in their study is the "median" SFR in a given 
mass bin, which is differe nt from the s patiall y "a veraged" SFR 
of the nearby galaxies in lElbaz et aP (l2007h and ICooper et al.l 
(120081) ■ lPatelet al.1 (l2009h argued that the contradictory results 
may be caused by the different environment that they investigate 
and by the different sample selection. They observed the central 
region of galaxy clusters (that is not well resolved by the other 
studies). Moreover, they used a mass-limited sample of galaxies 
that is different from the luminosity-limited sample of galaxies 
used in other studies, and the latter can be affected by the contri- 
bution of low-mass, blue star-forming galaxies. 

In this paper, we study the environmental dependence of 
(U)LIRGs found in the recent SDSS data to improve under- 
standing of the triggering mechanism of (U)LIRGs' SFA and 
of the evolution of the SFR-environment relation. Section |2]de- 
scribes the data and the environmental indicators used in this 
study. Environmental dependence of the physical properties is 
given in |3] Discussions and conclusions are given in ^and ^ 
respectively. Throughout, we adopt h - 0.7 and a flat ACDM 
cosmology with density parameters Qa = 0.73 and Q,„ - 0.27. 



2. Data 

2.1. Observational Data Set 

To identify LIRGs, we first cross-correlate the largest available 
IR source catalog and redshift catalog. For the IR source catalog, 
we us e the IRAS Faint Sources Catalog - Version 2 (Moshir et al. 
119921 hereafter FSC92), which contains 173,044 IR sources at 
I b \> lOdeg having measured fluxes at 12, 25, 60 and 100 fim. 



For the redshift catalog, we use a spectroscopic sample 
of galaxies including the main galaxy sample (m^ < 17.77) 
and faint galaxies (m, > 17.77) wh ose spectroscopic red - 
shifts are available in the SDSS DR7 (lAbazaiian et al.ll2009l) . 
Completeness of the spectroscopic data in SDSS is poor for 
bright galaxies with nir < 14.5 because of the problems 
of saturation and cross-talk in the spectrograph, and for the 
galaxies located in high-density regions such as galaxy clus- 
ters due to the fiber collision. Thus, it is necessary to sup- 
plement the missing galaxy data to reduce the possible ef- 
fects of the incompleteness problem. Therefore, we also use 
a photometric sample of galaxies with m,. < 17.77 whose 
redshift information is available in the literature. We down- 
load the photometric sample of galaxies with nir < 17.77 and 
match it with the redshift catalogs of galaxies such as Updated 
Zwicky Catalog (UZC; Falco et al. 1999), MX northern Abell 
cluster re dshift survey ( Sli nglend et alj Il998li) . EFAR survey 
("Wegner et al. 1996. .19 99). the redshift survey of Abell clusters 
(Hill & Oegerlel ll993r [l998). Furthermore, we add redshift in- 
formation of the photometric sample of galaxies located within 
ten times the virial radius of some galax y clusters from the 
NASA Extragalactic Database (NED) (Park & Hwang 2009). In 
total, the redshift information for 184,845 galaxies in the photo- 
metric sample is compiled, which is overlapped with 173,828 
galaxies in the spectroscopic sample. Finally, we add 11,017 
galaxies to the spectroscopic sample of galaxies, which yields a 
final sample of 929,234 galaxies. In the result, the spectroscopic 
completeness of our sample is higher than 85% at all magnitudes 
with nir < 17.77 and even in t he center of galaxy clusters (see 
Fig. 1 in iPark & Hwan g"2009'). It is to be noted that the final 
sample also contains the galaxies whose spectral classifications 
provided in the SDSS, are QSOs. 

In order to investigate the physical parameters of galaxies, 
we use several value-added galaxy catalogs (VAGCs) drawn 
from SDSS data. Photomet ric and structure par ameters are 
adopted from SDSS pipeline (IStoughton et alJl2002h . The veloc- 
ity dispersion of gala xies is adopted from DR7 release of VAGC 
dBlanton et alj l20()5h . The stellar mass estimates are obtained 
from MPA/JHU DR7 VAGCQ, which is based on the fit of SDSS 
five-ba n d photometry with th e model o f Bruzua l & Charloj 
(I2003h dKauffmannetaLl 120031; lGallazzietal...2005 i The Ha 
fluxes are also taken from MPA/JHU DR7 VAGC, which 
is computed from Gaussian fit to continuum subtracted data 
(Tremontiet al. 2004). 

We adopt the galaxy morphology information from Korea 
Institu te for Adva nced Study (KIAS) DR7 VAGC ( Park & Choil 
2001; IChoi et al.l [2007; Choi 2010 in prep.), which contains 
697,320 main galaxy sample in NYU VAGCs as well as 10,497 
photometric sample of galaxies with the redshift information 
from various existing redshift catalogs. We perform additional 
visual classification for the galaxies in DR7 that are not included 
in KIAS DR7 VAGC. During the visual inspection of color im- 
ages of galaxies, we eliminate 403 spurious sources (e.g., faint 
fragments of bright galaxies, diffraction spikes of bright stars), 
and they are not included in the final sample of 929,234 galaxies. 

2.2. The Cross-Correlation 

The positional uncertainties of the IRAS sources are different de- 
pending on the scan direction (typically 5" in-scan direction and 
16" for the cross-scan direction, 1 cr), and vary from source to 
source (1"-13" for the in-scan direction, or minor axis of posi- 
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Fig. 1. Comparison of IR luminosities (Lir) computed from the 
IRAS 60 i-im fluxes only (LiR(from bO/im)) with those from the IRAS 
60 and 100 fim (LiR(fi-oin 60,ioo/im))- We only show 30% of the total 
sample. 



tional uncertainty ellipse and 3"-55" for the cross-scan direc- 
tion, or major axis of positional uncertainty ellipse, 1 cr). The 
positional uncertainties of the optical counterparts are negligible 
in comparison. If a galaxy in SDSS lies within 3 cr positional 
uncertainty ellipse of the IRAS source, we regard it a match. It 
gave us 15,611 IRAS sources having optical counterparts in the 
SDSS. When there are more than one SDSS galaxy within the 
IRAS ellipse (10% of the total number of matchings), we chose 
the one closest to the center of the IRAS source. 

We use the spectral energy distribution (SED) models of 
ICharv & Elbad (| 2001 ) to compute the IR luminosity of galaxies. 
We restrict the computation to 15,547 sources whose IRAS 60 
jum fluxes are reliable, which means that the flux quality flags are 
either 'high quality' or 'moderate quality'. Fluxes at 60 and 100 
fim are used for the SED fit when their fluxes are reliable. If 100 
fim flux is not reliable, we use only 60 fj.m flux for the SED fit. 
It was pointed out that 100 Atm flux of some IRAS sources could 
be overestimated (^ Jeong et al1l200 7'). but its effect on the estima- 
tion of IR luminosity is expected to be less than a factor of two. 
Since flux quality flags are different depending on the sources, 
we show how the derived IR luminosities agree with each other 
in Fig. [U It shows that IR luminosities computed from the IRAS 
60 and 100 fim (LiR(fi-om 60,i()0/im)) and those from the IRAS 60 yum 
fluxes only (LiR(fi-om 60//m)) agree well with rms=0.08 in a loga- 
rithmic scale. We further reject 40 galaxies whose redshift confi- 
dence parameters (zConf) are smaller than 0.65, and finally have 
a sample of 15,507 IR galaxies. 

Since our sample is constructed by the cross-correlation 
based on the position alone, it may contain spurious sources due 
to the chance presence of a galaxy within the positional uncer- 
tainty of IR sources. In order to secure the sample of highly re- 
liable association between IR sources and the optical counter- 
parts, we compute the likelih ood ratio (LR) for each association 
(ISuflierland & Saundersl[T992) . The LR is defined by the ratio of 
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Fig, 2, (a) Distribution of likehhood ratio (LR). (b) Reliability 
of our cross-correlation vs. LR. Vertical dashed line represents 
the critical LR values (90%) of reliable identification for IRAS 
sources adopted in this study. 



the probability of a true association to that of a chance associa- 
tioifl 



LR 



Q(< m) exp(-r^/2) 
27rcri(T2«(< m) 



(1) 



where n(< m) is the local surface density of objects brighter 
than the candidate. The "normaUzed distance" r is given by 



0-2 



(2) 



where di and are positional differences along the two axes of 
an error ellipse for an IR source, and cr[, 0-2 are the lengths of 
these axes. Since the positional uncertainties of SDSS galaxies 
are negligible compared to those of IRAS sources, we define cr 
as the length of the error axes of the IRAS sources. Q{< m) is 
a multiplicative factor which is a prior probability that a "true" 
optical counterpart brighter than the magnitude limit exists in the 
association, and for simplicity we set 2 = 1 in this study. 

We compute n(< m) using the photometric sample of galax- 
ies within 3 cr error ellipses. 



n(< m) — 



N(< m) 



(3) 



where A^(< m) represents the number of galaxies of which mag- 
nitudes are less than or equal to that of a candidate, and a = 3. 
We then obtain the LR for our sample 



LR = 



exp(-r/2) 
2A^(< m) 



(4) 



We compute LR values for all IRAS sources that have optical 
counterparts in SDSS using the photometric sample of galax- 
ies. In order to calculate the reliability of the association using 
the LR values, we perform random associations by offsetting the 
positions of IRAS sources by x 30', and recompute LR values 
for each random association (e.g.. iHwang et al..,2007 ). Using the 



We assume that the positional errors of IRAS sources are Gaussian. 
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distribution of LR values for true and random associations, the 
reliability of each association with a given LR is defined by 



10 100 
SFE„ (Mg yr-') 

Fig. 3. Comparison of SFRs determined from Lir with those 
from Luo- We only plot galaxies whose spectral type is classi- 
fied as star-forming (see 33.1b . Filled circles in (b) are median 
values of the ratios (SFRho/SFRir) at each bin, and the error 
bars indicate the dispersion of the ratios. 




a.OO 0,05 0,10 0.15 0.20 0,25 0,30 

Fig. 4. (a) IR luminosity, (b) stellar mass, and (c) absolute r- 
band magnitude vs. redshift. Red dots indicate IRGs (only 30% 
of IRGs are shown), and black dots denote galaxies without IR 
detection in the spectroscopic sample (only 1 % of galaxies in 
the total sample are shown). Solid and dashed lines in (c) define 
the volume limited samples, LI and L2, respectively. The bottom 
curve corresponds to the apparent magnitude limit (nir = 17.77) 
for the main galaxy sample in SDSS. 



RiLR) = 1 - 



„(LR) 



Mme(LR) ' 



(5) 



where Aftrue(LR) and Mandom(LR) are the number of true and ran- 
dom associations with a given LR. In Fig.|2] we present the dis- 
tribution of LR values and the reliability of our cross-correlation 
as a function of LR. It appears that the reliability for IRAS 
sources decreases as LR decreases. Therefore, we made a final 
sample of 13,470 IR galaxies by selecting IRAS sources whose 
reliability is larger than 90% (LR> 0.55). 

In the present study, we assign a single optical counterpart 
to a given IRAS source. However, there can be multiple op- 
tical sources within the 3 cr positional uncertainty ellipse of 
the IRAS source, which may contribute to the measured IR 
flux. This may result in the overestimation of Lir for the op- 
tical source, with larg er effect in high-density regions such as 
galaxy groups/cluste rs (lAUam et al.lll996l:lLe Floc'h et al.ll2002l : 
iBitsakiset all 12010 '). To investigate this effect on our galaxy 
sample, we checked the ratio of the IR luminosity (Lir) to the op- 
tical luminosity in r-band (L^) as a function of the background 
density and the distance to galaxy clusters for galaxies in the 
volume-limited sample of LI (to defined in the next section). 
IR luminosity is known to correlate with the optical luminos- 
ity (see Fig. |5^). If we had overestimated Lir for galaxies in 
high-density regions, the ratio Lir/L,- would be expected to be 
artificially larger in high-density regions. Instead, we find that 
this ratio does not change with the background density and the 
distance to galaxy clusters, which suggests that this effect is not 
significant here. The AKARJ satellite is going to provide new 
all-sky maps with a better spa tial resolution than IR AS, which 
will definitely solve this issue (Murak ami et al.ll2007h . 

From now on, we call galaxies with Lir > IO'^Lq infrared 
galaxies (IRGs), those with 10'° < Lir < 10'°^Lo lower- 
luminosity infrared galaxies (LLIRGs), those with 10'*'^ < 
Lir < IO'^Lq moderately luminous infrared galaxies (MLIRGs), 
those with 10" < Lir < IO'^Lq LIRGs, those with 10'^ < Lir < 
IO'^Lq ULIRGs, and those with Lir > lO'^L© hyperluminous 
infrared galaxies (HLIRGs) (Goto 2005). In the final sample, we 
have 11, 964 IRGs including 3049 LLIRGs, 4362 MLIRGs, 4266 
LIRGs, 268 ULIRGs, and 19 HLIRGs. 

Since the IRGs in this study have IR luminosities and Ha 
emission line fluxes provided by SDSS, we compared the two 
SFR measurements for further analysis . We convert the I R lumi- 
nosity into SFRiR using the relation in lKennicut t'(n^9?): SFRir 
(Mq yr"') = 1.72 X 1O""'Lir(L0). The Ha luminosity (Lna) is 
computed from the flux measurement given in the MPA/JHU 
VAGCs by applying aperture corrections following the method 
given in appendix A of Hopkins et al. (2003). Obs curation cor- 
rection is also applied using the reddening curve of lCardelli et al.l 
(Il989l) . Thus, SFR/fg is converted from the Lna using the rela- 
tion in iKennicuttI (119981) : SFRna (Mq yr ') = 7.9 x 10 ^^^^h^ 
(ergs s"'). Fig.[3]represents the comparison between SFRh^ and 
SFRiR, and shows that the two measurements seem to agree well 
although with a large r scatter as compared to other studies lim- 
ited to local galaxies d Kewlev et ani2002t iKennicutt et alJl2009h 
most probably due to the inclusion of more distant galaxies here. 
It is also seen that SFR//a^ tends to underestimate the SFR for IR 
luminous galaxies probably due to the increasing uncertainty on 
the obscuration correction, which is also seen in other studie s 
(iKewlev et al.ll2002l:lHopkins et al.ll2003HKennicutt et al.ll2009l) . 
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Fig. 5. Physical parameters of IRGs as a function of Lir. (a) ab- 
solute magnitude Mr, (b) (u - r) color, (c) stellar mass Mstar, 
and (d) velocity dispersion cr. Filled circles are median values of 
each bin, and error bars indicate the dispersion of the data, (e) 
AGN fraction (Seyferts, LINERs and composite galaxies). We 
show only 30% of IRGs in the plot, but median values and AGN 
fraction are computed using the whole IRG sample. 
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Fig. 6. (a) Lir of IRGs in LI as a function of E5. Filled circles 
are median values at each density bin. (b) Fraction of LIRGs 
plus ULIRGs among IRGs as a function of S5 . (c) Fraction of 
LIRGs plus ULIRGs among SF galaxies in LI (SFRh,, > 1 or 
SFRiR > 1 Moyr '). (d) Fraction of LIRGs plus ULIRGs among 
all galaxies in LI. (e) Fraction of SF galaxies among all galaxies 
in LI. 



2.3. Galaxy Environment 

We define a volume-limited sample of galaxies selected by ab- 
solute magnitude and redshift limits as shown in Fig. EJc). We 
define two samples, LI and L2, that were chosen to maximize 
the number of target IRGs (to be discussed in the end of this sec- 
tion) and the number of ULIRGs in the volume, respectively: LI 
(0.03 <z< 0.0788 and < -19.91) having 134,272 galaxies 
including 5513 IRGs (6 ULIRGs) and L2 (0.03 < z < 0.17 and 
Mr < -21.68) having 129,851 galaxies including 2157 IRGs (38 
ULIRGs). The apparent magnitude limit line (m^ - 17.77) for 
the SDSS main galaxy sample shown in Fig.|4fc) is o btained us- 
ing the mean /iT-correction relation given by eq. (2) of lChoi et alj 
(I2007h . 

We consider three kinds of environment indicators: a surface 
galaxy number density estimated from the five nearest neighbor 
galaxy (E5) as a large-scale environmental parameter, a distance 
to the nearest neighbor galaxy (/?„) as a small-scale environmen- 
tal parameter, and a distance to a galaxy cluster (R). 

The background density, S5, is defined by S5 - 5{nD^ ,) ^ 
Dp s is the projected distance to the 5th-nearest neighbor The 
5th-nearest neighbor of each target galaxy is identified in the 
volume-limited sample of LI among the neighbor galaxies that 
have velocities relative to the target galaxy less than 1500 km s"' 



to exclude foreground and background galaxies. The E5 for L2 
galaxies is computed using the same method as for LI galaxies. 

To investigate the environment of galaxy clusters in detail, 
we compute the distance to galaxy clusters as a second envi- 
ronme ntal parameter using the method given in iPark & Hwand 
(120091) . We use the Abell clusters and the sample of 929,234 
galaxies in this study. The radius of rioo^ci (usually called the 
virial radius) for each cluster where the mean overdensity drops 
to 200 times the critical d ensity of the un i verse Pc is computed 
using the formula given bv lCarlbergetan(ll997h : 

r2oo,ci = (6) 

where cTci is the velocity dispersion of a cluster and the Hubble 
parameter at z is H^(z) = H^[O.Jl + z)^ + Qi(l + zf + Qa] 
(Peebles 1993). Qm, Q^, and Qa are the dimensionless density 
parameters. 

To define the small-scale environmental parameter attributed 
to the nearest neighbor, we first find the nearest neighbor of a tar- 
get galaxy that is the closest to the target galaxy on the projected 
sky and satisfies the conditions of magnitude and relative veloc- 
ity. We searched for the nearest neighbor galaxy among galaxies 
that have magnitudes brighter than Mr = M^.target + 0.5 and have 
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Fig. 7. (a) Distributions of IRGs in LI in the space of Mstar and S5. ULIRGs, LIRGs, MLIRGs, LLIRGs, and non-SF galaxies are 
plotted with star, circle, triangle, upside down triangle, and square, respectively. Non-SF galaxies are those with SFRhq^ = and 
SFRiR - {Mq yr '). Error bars indicate the Icr-dispersion of the distribution. Grey contour map is the number density of non-SF 
galaxies, and red dots are ULIRGs. (b) /(U)lirgs contours in the space of Mstar and E5. Dots are IRGs. 
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Fig. 9. (a) LiR of IRGs in LI as a function of the clustercentric 
radius normalized to the cluster virial radius (^/r2oo,ci)- Filled 
circle is the median value of each density bin. (b) /(u)lirgs vs. 

■'^/'"200,cl- 



relative velocities less than Av - lyneighbois - ^targetl — 600 km 
s"' for early-type target galaxies and less than Av = 400 km s"' 
for late-type target galaxies. Thes e velocity l i mits cover most 
close neighbors as seen in Fig. 1 of iPark et al.1 (120081) . aince we 
use a volume-limited sample of galaxies with Mr < -19.91 (for 
LI) and < -21.68 (for L2), we often restrict our analysis 
to target galaxies brighter than M,-_target = -20.41 (for LI) and 
target — —22.18 (for L2) so that their neighbors are complete 
when we investigate the effects of the nearest neighbor galaxy. 



The virial radius of a galaxy within which the mean mass 
density is 200 times the critical density of the universe (p^) is 
computed by 

r„r = (3yL/4ny'\200p,)-''\ (7) 

where L is the galaxy luminosity, and y is the mass-to-light 
ratio. We assume that the mass-to-light ratio of early-type galax- 
ies is on average twice as large as that of late-type galaxies at the 
same absol ute magnitude Mr, w hich mea ns y(early)- 2y(l ate) 
[see §2.5 of iPark & Choi (l2009h and §2 of lParket al.i(l2008M . 
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Fig. 10. (a) Distributions of IRGs in LI in the space of the stellar mass and the clustercentric radius. Symbols are same as in Fig.|7] 
(b) /(U)LiRGs contours in the space of the stellar mass and the clustercentric radius. Dots are IRGs. 



Since we adopt f2,„ = 0.27, 200pc - 200p/D.m - 740p where 
p is the mean density of the universe. The value of mean density 
of the universe, p = (0.0223 + 0.0005)(yL)_2o(/^"'Mpc)-^ was 
adopted w here (yL)^20 is t he mass of a late-type galaxy with 
Mr = -20 (i Park et al.ll2008l) . According to our formula the virial 
radii of galaxies with Mr = -19.5, -20.0, and -20.5 are 260, 
300, and 350 /z"' kpc for early types, and 210, 240, and 280 /i"' 
kpc for late types, respectively. 



3. Results 

3.1. Physical Parameters as a function of IR luminosity 

We show several physical parameters as a function of IR lumi- 
nosity for IRGs in Fig.|5] In this plot, we use all the IRGs instead 
of the volume-limited sample of galaxies. It is seen that, as Lir 
increases, IRGs become more luminous in r-band (a), bluer in 
(m - r) color (b) and more massive (c, d). Interestingly, bluer 
colors for more IR luminous galaxies may be incompatible with 
the idea that IR luminous galaxies are dust-enshrouded systems. 
However, this can be explained by the idea that the optical colors 
of ULIRGs are dominated by the extended distri bution of young , 
blue stars, and the dust extinction is patchy (Chen et aPboiOh . 
[Veilleux et al. (1999) also found bluer optical colors for more IR 
luminous galaxies in the sample of ULIRGs, and suggested that 
the optical continuum was less affected by dust extinction than 
the emission-line gas near hot stars or AGNs. 

We measure the fraction of galaxies showing AGN activity 
among emission-line galaxies in (e). We determine the spec- 
tral U[ges_of£mission^line galaxies based on the criteria given 
by 'Kewley et al.' ('2006) using the emission line ratio diagram, 
commonly^own as Baldwin-Phillips-Terlevich (BPT) diagram 
dBaldwin et al.lll98li) : star-forming galaxies, Seyferts, LINERs, 
composite galaxies, and ambiguous galaxies. We regard Seyferts 
plus LINERs and composite galaxies as galaxies showing AGN 
activity, and restrict our analysis to the galaxies at z > 0.04 due 
to the problem of small (3") fixed-size aperture (iKewley et al.l 
|2006|) . The uncertainties of the fraction represent 68% (Icr) con- 
fidence intervals that are determined by the bootstrap resampling 



method. Interestingly, the AGN fraction for IRGs keeps increas- 
ing as LiR increases over all the IR luminosity range, which is 
consistent with p revious results (e.g.. lVeilleux et aP 1995Ll2002t 
lYuanet al.ll201(i^ . 

We wish to emphasize that the sample of (U)LIRGs used 
in this study is dominated by LIRGs (99.5% and 97.4% among 
(U)LIRGs in LI and L2, respectively), thus the results in next 
sections are mostly determined by LIRGs (to be discussed in 



3.2. Background Density Dependence of IRG Properties 

In Fig. |6] we present the IR luminosity of IRGs and the frac- 
tion of LIRGs plus ULIRGs among several galaxy samples 
(/(U)LiRGs) in LI as a function of the background density (S5). 

It appears that neither the median Lir values (Fig. |6h) nor 
/(U)LiRGs (in the IRG sample, b) change with Z5. Note that 
/(U)LiRGs in (b) is the fraction of LIRGs plus ULIRGs among 
IRGs, which is dififerent from the fraction of SF galaxies that 
are usually found to decrease monotonously as the background 
density increases. This is due to the large contribution of non- 
SF galaxies in high-density regions (see e). When we compute 
/(U)LiRGs in the whole sample regardless of IR and Ha detections 
(Fig.|6}l), it decreases as 25 increases. On the other hand, if we 
compute /(U)LiRGs in the SF galaxy sample as shown in (c), we 
find a similar trend of /(u)lirgs to the case of IRGs in (b). We 
call SF galaxies those with SFRir or SFRh^ > 1 Moyr We 
include galaxies having Ha measurements in order to take into 
account for SF galaxies that are not detected in IR bands due to 
the IR detection limit. Since IRGs can be regarded as a repre- 
sentative sample of SF galaxies (IRGs are by definition galaxies 
with SFRiR > 1.72 Moyr"' in this study) and we are interested 
in the triggering mechanism of (U)LIRGs that are extreme cases 
among IRGs, we are going to focus on /(u)lirgs in the IRG sam- 
ple in further analysis. When we use only the sample of IRGs 
hosting AGNs in (b), no dependence of /(U)lirgs on the back- 
ground density is seen again. 
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Fig. 11. (a) LiR of IRGs in LI plus L2 (-20.41 > M, > -23.76) 
as a function of the projected distance to the nearest neigh- 
bor galaxy (^n/ '"vir.nei)- Red and blue dots are the galaxies hav- 
ing early- and late-type nearest neighbor galaxies, respectively. 
Circle (early-type neighbor case) and square (late-type neighbor 
case) are median values of each density bin. /(u)lirgs vs. rvi,- n 
(b) for IRGs in LI plus L2 (-20.41 > M, > -23.76) and (c) for 
IRGs with -21.20 > > -21.80 among LI plus L2. 



Since the SFA of galaxies is s trongly linked to th e 
stellar mass (Bri nchmann et alj l2004t iNoeske et alj l2007bllal: 
lElbaz et a l. 2007), it is important to quantify its effect in or- 
der to investigate the effects of environment. Thus we plot the 
distribution of IRGs in company with the non-SF galaxies in 
two-dimensional space of the stellar mass (Mstar) and the back- 
ground density (£5) in Fig.|3a). Non-SF galaxies are those hav- 
ing neither Ha line emissions nor IR emission (SFRhu^ - and 
SFRiR - Mq yr"'). It is clearly seen that the stellar mass of 
IRGs increases as Ljr increase, while the median values of do 
not change much with Lir. IRGs tend to be found in low-density 
regions compared to non-SF galaxies. Panel (b) shows constant 
/(U)LiRGs contours in the same parameter space. Interestingly, the 
/(U)LiRGs contours are nearly horizontal, which confirms that the 
dependence of /(u)lirgs on the background density is very weak 
compared to that on the stellar mass. When we use galaxies in L2 
that contain more ULIRGs than LI, no dependence of /(u)lirgs 
is similarly found again (see Fig. [8]). 

3.3. Effect of Galaxy Clusters on IRG Properties 

In Fig. |9] we present Lir of IRGs and /(u)lirgs as a function 
of the projected distance to the galaxy cluster {R). We restrict 
this plot to the IRGs whose clustercentric distances are mea- 
sured in this study. Fig. |9ja) shows that the median Lir values 
do not change much with the clustercentric distance. In panel 



(b), /(u)LiRGs appears to decrease with large error-bars, as we ap- 
proach the center of galaxy clusters. Moreover, in the very center 
of galaxy cluster (R < 0.5r2oo,ci), no (U)LIRGs are found. 

We plot, in Fig. (TO] the distribution of IRGs in the space of 
the stellar mass and the clustercentric radius. The right panel 
shows that /(u)lirgs decreases as IRGs are located closer to the 
center of clusters. The dependence of /(u)lirgs on the cluster- 
centric radius is much weaker than the dependence on the stellar 
mass. According to our finding in the next section this depen- 
dence can appear because the typical morphology of the nearest 
neighbor changes from late type to early type. The IRGs near the 
central region are typically under the strong influence of early- 
type neighbors, and are diflicult to show (U)LIRGs behaviors. 

3.4. Nearest Neighbor Dependence of IRG Properties 

To measure the effect of the nearest neighbor galaxy on the SFA 
of IRGs, we plot Lir of IRGs and /(u)lirgs as a function of the 
distance to the nearest neighbor in Fig. [TT] In order to remain 
complete in terms of neighbor galaxies, we restrict our analysis 
to the galaxies having M,- < -20.41 and Mr < -22.18 in LI 
and L2, respectively. Panels (a-b) show that Lir and /(u)lirgs de- 
pend on the projected distance to the nearest neighbor as well 
as on the neighbor's morphology. When an IRG is located far- 
ther than 0.5 rvir,nei (virial radius of neighbor galaxy), /(u)lirgs 
decreases as the distance to the neighbor decreases, but its de- 
pendence on the neighbor's morphology is negligible. On the 
other hand, when an IRG is located at Rn < 0.5rvii,nei, /(U)lirgs 
increases as the target IRG approaches a late-type neighbor, but 
decreases as it approaches an early-type neighbor It is important 
to note that the bifurcation of /(U)lirgs depending on the neigh- 
bor's morphology occurs at Rn ~ 0.5rvii,nei, which is similarly 
seen in the plot of median Lir values along with the distance to 
the nearest neighbor. 

In order to check whether the difference in /(u)lirgs between 
IRGs with a late- and early-type close neighbor is real or if it 
could be explained by small number statistics, we made the fol- 
lowing experiment. We produced 1000 Monte-Carlo realizations 
of the same number of galaxies with close neighbors by ran- 
domly selecting galaxies from the larger sample of galaxies with 
no neighbor closer than Rn = '"vir.nei (24 and 7 galaxies with late- 
and early-type neighbors at R„ < 0.03rvii,nei, respectively). While 
in the real dataset, /(U)lirgs is 3.8 times larger for IRGs with a 
late-type neighbor than for those with an early-type neighbor, 
such large difference is found in only 6% of the random samples 
which suggests that the difference is real. This experiment sup- 
ports the robustness of the bifurcation of /(u)lirgs depending on 
the neighbor's morphology at R„ < O.Sryir.nei- 

The trends seen in panel (b) of Fig. [TT]with respect to R„ 
are actually coupled with luminosity dependence of /(u)lirgs ■ To 
disentangle /(u)lirgs from the coupling with luminosity and to 
see its dependence on the neighbor environment only, we fix the 
absolute magnitude of IRGs to a narrow range of -2 1 .20 > M, > 
-21.80 in panel (c) of Fig.[TT] Now /(u)lirgs remains nearly con- 
stant down to Rn ~ 0.5rvir,nei, which indicates no direct neighbor 
effect at these separations. The bifurcation of /(U)lirgs still oc- 
curs at Rn ~ O.Srvir.nei- We will discuss these findings in more 
detail in S4J\ 

In Fig. [12] we present the distribution of IRGs in the space of 
the stellar mass and the distance to the nearest neighbor When 
the nearest neighbor galaxy of IRGs is an early type (a), IRGs 
are preferentially located farther away from their neighbor than 
non-SF galaxies, and more luminous IRGs are farther away from 
their neighbor. On the other hand, the median distance to the 
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Fig. 12. Distributions of IRGs in LI plus L2 in the space of the stellar mass and the distance to the nearest neighbor galaxy when 
the nearest galaxy is an early type (a) and a late type (b). Symbols are same as in Fig.]?] Red dots are ULIRGs. 
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Fig. 13. /(u)LiRG.s for IRGs in in LI plus L2 in the space of Mstar and 7?n/'"vir,nei when the nearest galaxy is an early type (a) and a late 
type (b). Dots are IRGs. 



nearest neighbor for the late-type neighbor case (b) is not clearly 
distinguishable among IRGs including non-SF galaxies. 

The /(U)LiRGs contours in the same parameter space are plot- 
ted in Fig. [T3] depending on the morphology of the nearest 
neighbor. Unlike the /(u)lirgs contours in the space of Mstm 
and ^5 seen in Fig. the constant /(u)lirgs contours in Fig. 
[T3] depend on both M^tar and When an IRG is located at 
Rn ^ 0.5rvir,nei, the probability for IRGs to become (U)LIRGs in- 
creases as IRGs are more massive and are farther away from their 
neighbor regardless of the morphology of the nearest neighbor. 
However, if an IRG is located at ^ 0.5rvir,nei, the dependence 
of /(u)LiRGs on Mstar and R„ is significantly different depending 
on the morphology of the nearest neighbor. For the early-type 
neighbor case, /(U)lirgs increases as IRGs are more massive and 



are farther away from their neighbor, while /(U)lirgs increases as 
IRGs are more massive and are closer to their neighbor for the 
late-type neighbor case. The different behavior of /(U)lirgs de- 
pending on the morphology of the nearest neighbor indicates an 
important role of the nearest neighbor on the activity of IRGs. 

Bearing in mind that there is a strong effect of the nearest 
neighbor galaxy on /(U)lirgs, we re-investigate the dependence 
of /(U)LiRGs on the background density and the distance to galaxy 
clusters. Fig.[T4]shows /(U)lirgs as a function of the background 
density and the distance to galaxy clusters as seen in Figs. |6] 
|8] and |9] but depending on the morphology of nearest neighbor 
galaxy. It is seen that /(u)lirgs for the late-type neighbor case, on 
average, are larger than those for early-type neighbor case. The 
difference of /(U)lirgs are insignificant in low-density regions 
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Fig. 14. (a) /(u)LiRGs for IRGs in LI as a function of E5. Circle 
(early-type neighbor case) and square (late-type neighbor case) 
are median values of each density bin. Vertical dashed line indi- 
cates the Z5 that roughly corresponds to one virial radius of the 
clustercentric radius (R ~ r2oo,ci) determined from the compar- 
ison between E5 and the distance to galaxy clusters in LI. (b) 
Same as above, but for /(u)lirgs vs. /J/raoo.ci- 



where the distance to the nearest neighbor is large. However, 
the difference of /(u)lirgs becomes prominent in high-density 
regions (even in the cluster region of R ~ 0.6r2oo,ci) where 
the neighbor separation is smaller compared to low-density re- 
gions so that IRGs experience the effect of the neighbor galaxy. 
These confirm again the strong effects of the nearest neighbor on 

/(U)LIRGs- 

4. Discussions 

4.1. What powers LIRGs and ULIRGs? 

The strong dependence of /(u)lirg,s and median Lir values on 
the morphology of and the distance to the nearest neighbor 
galaxy seen in ^3.41 supports the idea that LIRGs and ULIRGs 
are triggered by galaxy-galaxy interactions/merging. In particu- 
lar, the bifurcation of /(u)lirgs and median Lir values occurs at 
Rn ~ 0.5rvir,nei where th e effects of galaxy-galaxy interactions 
can significantly appear (iPark & Choi 2009) . The bifurcation of 
/(U)LiRGs and median Lir values at Rn ~ 0.5rvir,nei depending on 
the morphology of the neighbor may imply that the hydrody- 
namic interactions with the nearest neighbor play a critical role 
in triggering SFA of IRGs in addition to the tidal interactions. 
If an IRG approaches a late-type neighbor within half of the 
virial radius of the neighbor, the inflow of cold gas from the 
neighbor into the target IRG increases and the SFA of IRGs can 
be boosted. On the other hand, if an IRG approaches an early- 
type neighbor within the virial radius of the neighbor, the hot gas 
of the early-type neighbor prevents the IRG from forming stars 
with cold gas or there is no inflow of cold gas from its early- 
type neighbor, so the SFA of IRG is not boosted even if it has 
a close companion. The SF quenching mechanisms of hot gas 
in early-type neighbors could be similar to those of hot intra- 
cluster medium of galaxy clusters acting on late-type galaxies in 
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Fig. 15. (a) Stellar mass ratio between IRGs and their nearest 
neighbor galaxy, and (b) projected distance to the nearest neigh- 
bor galaxy of IRGs (/?n/'"vir,nei) as a function of Lir of IRGs in 
LI plus L2. IRGs {Mr < -21) having late-type neighbors at 
Rn < '"vii.nei are shown. Large filled circles and vertical error- 
bars represent the median values and their errors in Lir bins that 
are represented by horizontal error-bars . The dashed horizontal 
line in (a) divides IRGs into major (mass ratio< 4) and minor 
(mass ratio> 4) mergers. 



it, which are hydrodynamic processes such as thermal evapora- 
tion, strangulation, ram pressure stripping, or viscous stripping 
dPark & ChoiH2009l:lPark & H wang 2009). IRGs may not neces- 
sarily be late types. In fact, 4% of IRGs in this study are early 
types, which shows again the importance of hydrodynamic in- 
teractions with late-type neighbors for the SFA. 

Interestingly, if we see Fig.[TT](a-b), /(u)lirgs and median Lir 
values at the largest neighbor separations (/?„ ~ 15rvii-,nei) appear 
to be larger than those at the intermediate neighbor separation 
{Rn ~ 2rvir,nei), which may provide another hint for the evolu- 
tion of (U)LIRGs. As an IRG experiences interactions/merging 
with its closest neighbor galaxy, the galaxies in pair tend to be 
(U)LIRGs though the probability to become (U)LIRGs is dif- 
ferent depending on the morphology of the neighbor When two 
galaxies just finished to merge, the end product of the merger 
will be bright due to the very recent SFA and/or the increase 
of mass, and the new nearest neighbor galaxy of the merger 
product will be far away. This seems to be related to the fact 
that some (U)LI RGs ar e in pair and othe rs are found as a sin- 
gle system (e.g.. lKim et al. 2002; Veilleux "etan l2002h . and to 
the variation of SFA of (U )LIRGs along with the merger stage 
(e.g., Murphv et al.ll2001h . The SFA then would decrease until 
the galaxy experiences new interactions/merging with the near- 
est neighbor galaxy. This scenario is consistent with the fact that 
there is no dependence of /(u)lirgs and median Lir values on 
the neighbor's morphology at this separation (/?„ > 0.5rvii,nei), 
and could explain why /(u)lirgs and median Lir values at the 
largest neighbor separation are larger than those at the interme- 
diate neighbor separation. As we show in Fig. fTTT c). if we fix the 
absolute magnitude of IRGs (or stellar mass), the bifurcation of 
/(U)LiRGs depending on the morphology of the nearest neighbor at 
Rn < 0.5rvir,nei IS Still sccn, but /(u)LiRGs (and median Lir values) 
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Fig. 16. (a) SFRs of IRGs in LI, (b) stellar mass of IRGs, Mst ai, and (c) specific SFRs of IRGs as a function of Sj. (d-f) Same as 
(a-c), but for the high-redshift IRGs (0.8 < z < 1.2) studied in lEIbaz et alJ (l2007h . Median value and its error at each bin is shown 
by filled circle and its error-bar. 



remains constant at large separations > 0.5rvii-,nei), which is 
consistent with the interpretation above. These results seem to 
be compatible with recent studies focusing on the effects of the 
nearest neighbor galaxies on the evolution of "normal" galajc- 
ies (Park et al. 2008; Park & Choii I2009t iPark & H wan"g 2009; 
iHwang & Park 2009) . and with the scenario of transformation of 
galaxy morphology a nd luminosity cla ss through galaxy-galaxy 
interactions/merging ( Park et al]|2008h . 

On the other hand, we found no dependence of /(u)lirgs on 
the background density (see Figs. |6l[8]l. When we use the dis- 
tance to the nearest galaxy cluster instead of using ^5, the re- 
sults are similar: dependence of /(u)lirgs on the clustercentric 
distance is weak except for highest-density regions such as the 
center of galaxy clusters. If the SFA of IRGs is increased by in- 
teractions in general, one would expect to find an increase of 
/(U)LiRGs with the background density. The fact that such trend is 
not seen is explained by the finding of a selective effect of the 
neighboring galaxy when it is closer than the virial radius of the 
neighbor, i.e. only late types produce an increase of /(u)lirg.s- 
This can be understood by the need to provide cold gas to- 
gether with the dynamical interaction. These results are compat- 
ible with the result from the study of Hickson compact groups 
of galaxies that specific SFRs of galaxies in dynamically 'young 
(dominated by late-type galaxies)' groups are, on average, larger 
than those in dynamically 'old (d ominated by early-type galax- 
ies)' groups (Bitsakis et al.'2010). 

As we mentioned in ^3.11 the sample of (U)LIRGs used in 
this study is dominated by LIRGs. Thus these findings demon- 
strate that not only ULIRGs but also LIRGs are strongly affected 
by interactions/merging. In order to investigate what makes 
LIRGs and ULIRGs different in detail, we plot the ratio of the 
stellar mass between IRGs (Mirg) and their nearest neighbor 
galaxy (Mn) (i.e. max[MiRG,M„]/mm[MiRG,M„]) as a function of 
LiR in Fig.fTSh. We use all galaxies in LI or L2 to find the nearest 



neighbor galaxy of IRGs without using the magnitude condition 
(Mr - Mr.target + 0.5) introduced in g2.3l in order to have no con- 
straints on the mass range of neighbor galaxies. We show only 
IRGs with Mr < -21 in order to reduce the effect of absolute 
magnitude limit in LI (Mr=-19.91) on the lower limit of neigh- 
bor's mass. In addition, we plot only IRGs interacting with a 
late-type nearest neighbor (i.e. IRGs with a late-type neighbor at 
^ '"vir.nei)- It shows that the mass ratio between an IRG and its 
interacting pair appears to decrease as Ljr increases. The fraction 
of minor merger system (i.e. mass ratio> 4) is found to be 36, 
34, and 14% for LLIRGs plus MLIRGs, LIRGs and ULIRGs, 
respectively. Due to the small number of ULIRGs in our sample, 
the fraction of minor merger system for ULIRGs is computed by 
including galaxies with Lir > 9 x 10"Lq. We check the statisti- 
cal significance of the difference in the fraction of minor merger 
system depending on Ljr by Monte Carlo simulations. We con- 
struct 1000 trial data sets by randomly selecting the same num- 
ber of galaxies at each Lir bin from the real data, and compute 
the fraction of simulated data sets in which the median value of 
mass ratio is larger than 4 (i.e. minor merger system). We find 17, 
18, and 2% for LLIRGs plus MLIRGs, LIRGs and ULIRGs, re- 
spectively. These results indicate that non-negligible fraction of 
LIRGs is related to the minor merger/interactions between two 
galaxies with a large mass difference, while few ULI RGs are mi- 
nor rn erger sy stems ((ish ida 2004; Wang et al. 2006; Dasvra et al 



2006; see also'Shi et al. 2006ti2009i;iElbaz et al.i2007l;lLotz et al 
2008; Ideue et al. 2009). 

Using the same sample as above, we show projected dis- 
tances to the nearest neighbor galaxy of IRGs as a function of 
Lir in Fig. [T5b . It is seen that distances to the nearest neighbor 
galaxies are smaller for ULIRGs than those for LIRGs. If IRGs 
with Rn < 0.1rvir,nei Can be regarded as advanced merger sys- 
tems, the fraction of advanced merger systems for LLIRGs plus 
MLIRGs, LIRGs and ULIRGs is 6, 21, and 50%, respectively. 
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Fig. 17. (a) Average SFR computed using all galaxies in LI, < aSFR >, (b) average SFR computed using IRGs, < SFR >, and 
(c) the fra ction of IRGs amo ng all galaxies as a function of E5. (d-f) Same as (a-c), but for the high-redshift IRGs (0.8 < z < 1.2) 
studied in lElbaz et al. I (l2007h . Median value and its error at each bin is shown by filled circle and its error-bar. 



We check the statistical significance of the difference in median 
values of ^n/'"vir,nei depending on Lir based on the Monte Carlo 
simulations with 1000 trial data sets made by randomly selecting 
the same number of galaxies at each Lir bin from the real data. 
We compute the fraction of simulated data sets in which the me- 
dian value of /?n/'"vir,nei IS Icss than 0.1 (i.e. advanced merger), 
and find 0, 3, and 49% for LLIRGs plus MLIRGs, LIRGs and 
ULIRGs, respectively. These results can suggest that ULIRGs 
are more advanced merger systems (or strong interaction) than 
LIRGs. 

In summary, the SFA of both LIRGs and ULIRGs is trig- 
gered by the interaction/merging with late-type (gas-rich) neigh- 
bor galaxies, but they seem to differ in the sense that 1) stellar 
masses of ULIRGs are larger than those of LIRGs (see Figs.|5}; 
and [121). 2) few ULIRGs are minor merger systems, but a non- 
negligible fraction of LIRGs (~34%) is a minor merger system, 
and 3) ULIRGs are more advanced merger systems than LIRGs. 



4.2. Evolution of SFR-environment Relation 

To investigate how the effects of the environment on the SFA of 
galaxies evolve with cosmic time, we plot, in Figs. [T6land[T7l 
several SFA parameters of IRGs as a function of the background 
d ensity (£5) in com parison with those at high redshifts studied 
in lElbaz etai] 120071 

In Fig. \W[ SFRs, stellar masses, and specific SFRs for local 
IRGs (a-c) do not show any significant change along with the 
background density, while SFRs and stellar masses for high red- 
shift IRGs (d-e) increase as the background density increases. 

In Fig.[T7] we show the spatially averaged SFRs for IRGs as 
a function of the background density. The average SFR is com- 
puted at the position of IRG using the neighbor galaxies that 
have velocities relative to the IRG less than 1500 km s"' and the 



projected distance to the IRG less than 1 Mpc. < aSFR > in (a 
and d) is the average SFR computed using all neighbor galaxies 
regardless of IR detection, while < SFR > in (b and e) is the one 
computed using the neighb or IRGs. For the GOODS galaxies 
studied in lElbaz et aTl ([2007), we first compute the background 
density (Sg.cooDs) by identifying 5th-nearest neighbor of each 
IRG in the spectroscopic sa mple of gala xies with 0.8 < z < 1.2 
and zab ^ 23.5 (see E lbaz et al.l 12007 * for details of the data). 
Secondly, we compute the average SFRs adopting the same 
method as above. 

As seen in (a) and (d), the average SFR (< aSFR >) for lo- 
cal SDSS IRGs decreases as the background density increases, 
while that for high redshift IRGs increases as the background 
density increases. This is the re versal of the SF-den sity relation 
in the distant universe found by lElbaz et al ] (I2007h . The differ- 
ent dependence of the average SFR on the background density 
might be caused by the different fraction of star-forming galax- 
ies depending on the background density (see c and f): since 
the fraction of star-forming galaxies in the local universe is very 
small in high-density regions, the average SFRs in these regions 
are also small due to the large contribution of non-SF galaxies. 

If we see the behavior of the average SFR (< SFR >) com- 
puted using the neighbor IRGs, < SFR > for local IRGs in (b) 
appears not to depend on the background density, while < SFR > 
for high redshift galaxies in (e) still shows an increase along 
with the background density. Since the average SFR (< SFR >) 
computed using the neighbor IRGs is not affected by the non-SF 
galaxies, it may provide important hints for what directly affects 
the SFA. Since the SFR of a galaxy shows a strong correlation 
with t he galaxy mass as seen in ^ (see also Brinchmann et al.l 
l2004t lElbaz et al.ll2007t iNoeske et al.ll2007 a"?r it is instructive 
to see the mass dependence of IRGs on the background density. 
It is seen in Fig. [TSl b) and (e) that the median stellar masses 
for local IRGs do not show a significant change with the back- 
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ground density, while those for high redshift galaxies increase as 
the background density increases, which are responsible for the 
difference of the average SFRs (< SFR >) between Fig. [TTI'b) 
and (e). 

Finally, the key in the SF-density relation seems to be how 
the fraction of massive, star-forming galaxies changes with the 
background density. Then the reversal of SF-density relation can 
be understood by the idea that massive, sta r-forming galaxies 
are strongly cl ustered at high redshifts (e.g., iFarrah et al.ll2006t 
iGilh et ani2()0 7) as forming a larger structure such as a galaxy 
group/cluster, and evolve (or consume their gas) through galaxy- 
galaxy interactions/merging faster than those in the field. Thus, 
gas-exhausted, non-SF galaxies are preferentially found in high- 
density regions of the local universe. 

This environmental dependence of evolution of the SFA 
of galaxies is consistent with the results of recent stud- 
ies o f other galaxy p roperties such as galaxy morphology 
(e.g.. ICapak etal.1 12007: Hwang & ParkI (2009) . For example, 
iHwang&ParkI (120091) found that the early-type fraction evolves 
much faster in high-density regions than in low-density regions, 
and that the morphology-density relation becomes weaker at 
Z ~ I (see their Fig. 8). It may be because the rate of galaxy- 
galaxy interactions is higher in high-density regions, and a series 
of interactions/merging over the course of galaxy life consume 
most of its cold gas and/or transform late types into early types. 

5. Conclusions 

By cross-correlating the spectroscopic sample of galaxies in 
the SDSS DR7 and the IRAS FSC92, we obtained a sample 
of nearby LIRGs and studied their environmental dependence. 
We examined the effects of the large-scale background density, 
galaxy clusters, and the nearest neighbor galaxy on the proper- 
ties of IRGs. Our main results are as follows. 

1. We find no dependence of the fraction of LIRGs plus 
ULIRGs among IRGs (/(u)lirgs) on the background density 
(S5), but find a strong dependence on stellar mass. 

2. When we use the distance to the nearest galaxy cluster in- 
stead of using ^5, the results are similar to those based on 
Z5: the dependence of /(U)lirgs on the distance to the nearest 
galaxy cluster is weak except for the highest-density region 
such as the center of galaxy clusters. 

3. We find that /(u)lirgs and Lir of IRGs strongly depend on 
the morphology of and the distance to the nearest neighbor 
galaxy. When an IRG has a close neighbor galaxy, the prob- 
ability for it to be a (U)LIRG (/(u)lirgs) and its Lir increase 
as it approaches a late-type neighbor, but decrease as it ap- 
proaches an early-type neighbor The bifurcations of Lir and 
/(U)LiRGs depending on the neighbor's morphology are seen 

at R„ ~ 0.5r vii-,nei- 

4. The median Lir values of IRGs and /(U)lirgs at the largest 
neighbor separations seem to be larger than those at the in- 
termediate neighbor separation. This supports the idea that 
(U)LIRGs are powered by the interaction/merging with late- 
type (gas-rich) nearest neighbor galaxies as expected for 
non-resolved advanced mergers. 

5. LIRGs and ULIRGs are found to differ in the sense that 1) 
stellar masses of ULIRGs are larger than those of LIRGs, 2) 
few ULIRGs are minor merger systems, but a non-negligible 
fraction of LIRGs (~34%) is a minor merger system, and 3) 
ULIRGs are more advanced merger systems than LIRGs. 

6. The environmental dependence of (U)LIRGs and the evolu- 
tion of the SFR-environment relation from high redshifts to 



low redshifts seem to be consistent with the idea that galaxy- 
galaxy interactions/merging play a critical role in triggering 
the SFA of (U)LIRGs. Massive, star-forming galaxies at high 
redshifts are strongly clustered, and evolve (or consume their 
gas) through galaxy-galaxy interactions/merging faster than 
those in the field. Thus, gas-exhausted, non-SF galaxies are 
preferentially found in high-density regions of the local uni- 
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